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Abstract 

We present a general relation between Mandel's Q param- 
eter describing variance of number of photons emitted from a 
single molecule and a three time correlation function describ- 
ing a spectral diffusion process the molecule is undergoing. 
For a prototype of such a stochastic process an exact solution 
is found which exhibits rich types of physical behaviors. 

Pacs: 78., 42.50.Ar, 33.80.-b, 05.40.-a 

Experimental advances have made it possible to mea- 
sure the spectral line of a single molecule (SM) embedded 
in a condensed phase Q . The statistical properties of the 
photons emitted from such a SM illuminates the interplay 
between dynamical processes in the condensed phase, SM 
dynamics and light-matter interaction. An important 
process responsible for the emission pattern of a SM is 
spectral diffusion, i.e., perturbers such as two level sys- 
tems or excitations in the environment of the excited SM 
lead to random changes in the transition frequency of the 
SM and in this way the SM comes in and out of resonance 
with the exciting laser field The fundamental is- 

sue addressed in this work is the photon statistics of a 
single molecule undergoing a spectral diffusion process. 
Photon statistics of an ensemble of molecules will clearly 
differ from that of a SM due to inhomogeneous broaden- 
ing which is absent in SM spectroscopy (SMS) and due 
to correlation between fluorescence photons which exists 
only on the SM level. 

We will model these fluctuations in a semiclassical way, 
using the Bloch equation, in the limit of weak laser field. 
The spectral diffusion process is described using a Kubo- 
Anderson sudden jump approach |12|. The same model 
leads to a basic understanding of line shape phenomena, 
namely of the average number of counts (n) per measure- 
ment time T. Variants of this model have been used in 
the past fifty years to model ensemble measurements in 
e.g., NMR Q and non-linear spectroscopy M] . More re- 
cently, this approach was applied to model SMS in glass, 
by Geva and Skinner jl4|] and in jllj] to describe static 
properties of line shapes and by Plakhotnik jlO| to model 
time dependent fluctuations of SMS. 

Here we present a general formula for Mandel's Q pa- 
rameter, 



Q 



1 



where n is the random number of counts. We show 
how Q is related to the underlying stochastic events and 



how Q yields new information not contained in the aver- 
age (n). For a simple Kubo- Anderson spectral diffusion 
process we find an exact solution for Q which gives in- 
sight on characteristic behaviors of SMS fluctuations. If 
Q = the photon count is Poissonian while our semi- 
classical results show a super Poissonian behavior Q > 0. 
For short times (e.g., ns), the correlation function be- 
tween fluorescence photons emitted by a single molecule 
shows anti-bunching, a sub-Poissonian non-classical ef- 
fect Our approach works well for longer times 
such that TT 3> 1 where T is the radiative decay rate of 
the molecule. 

The key quantity in the present formulation is a three 
time correlation function C3 (ti,T2,T3), which is similar 
to non-linear response functions investigated in the con- 
text of four wave mixing processes fl3|| . The three time 
correlation function contains all the microscopic infor- 
mation relevant for the calculation of Q. It appeared 
previously in a recent work of Plakhotnik [ic| l in the con- 
text of intensity-time-frequency-correlation technique. 
In this letter important time ordering properties of this 
function are investigated, and Q is evaluated. The re- 
lation between C3 {t\,T2,tz) and line shape fluctuations 
described by Q generalizes the Wiener-Khintchine theo- 
rem, which gives the relation between the averaged line 
shape and the one time dipole correlation function. 

An important issue is the fast and slow modulation 
limits, to be defined later. If the bath is fast the line 
exhibits the well known behavior of motional narrowing 
(as the time scale of the bath becomes short; the line is 
narrowed). Considering a simple spectral diffusion pro- 
cess, we show that Q also exhibits motional narrowing, 
though in this limit the fluctuations described by Q be- 
come small and when time scale of the bath goes to zero 
we find Poissonian statistics. If the bath is slow, Q be- 
comes large, and a simple adiabatic type of approxima- 
tion, based on steady state solution of time independent 
Bloch equation, holds. Our exact results are evaluated 
for an arbitrary time scale of the bath and are shown to 
interpolate between the fast and slow modulation lim- 
its. We note that several authors have analyzed intensity 
fluctuations in SMS using the slow modulation assump- 
tion , here we extend these results to include the 
physically important case of fast modulation. 

We assume a simple non-degenerate two state single 
molecule in an external classical laser field. The first 
electronic excited state |e) is located at energy ujq above 
the ground state \g). We consider the SM Hamiltonian 
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H = -u) (T z +'^2-Au j (t)<T z -d-'Ecos(u> h t), (1) 
i=i 

where o~ z is the Pauli matrix. The second term in H 
reflects the effect of the time evolution of the environ- 
ment on the absorption frequency of the SM coupled to 
J perturbers. The last term in Eq. (Q) describes the in- 
teraction between the SM and the laser field (frequency 
wl), while d = d cg a x is the dipole operator with the real 
matrix element d cg = (e|d|<?). 

The molecule is described by 2 x 2 density ma- 
trix p whose elements are p cc , p gg , p eg , p ge . Let, u = 
(p gc e-^t + /9eg e^*)/2, v = -i(Aj.e-*^* - p og e^*)/2, 
and w = {p cc — p gg )/2. Using Eq. (|l|) the stochastic Bloch 
equations, in the rotating wave approximation, are 

u = 6 h (t)v - Tu/2, 

v = —S^(t)u — flw — r^/2, 

w = Qv-Tw-T/2, (2) 

Q = — d cg ■ ~E/h is the Rabi frequency, 5^(t) = cu-^ — lu — 
Auij(t) is the detuning frequency and without loss 
of generality we set ujq = 0. SM experiments jl7|,[l8| were 
shown to be compatible with the deterministic two level 
optical Bloch equation approach and this gives further 
justification for our assumptions. 

According to semiclassical theory of photon counting 
statistics JL9|, the probability of recording n photons in 
time interval (0, T) is given by 



p(n,T) 



exp (-W) 



(3) 



with W = £ J Q I(t)dt, where I(t) is the photon current 
and £ is a suitable constant depending on detection effi- 
ciency. The semiclassical approach is valid for large num- 
ber of counts, when the stream of incoming absorption 
photons (i.e., Civ) is equal to the stream of emitted fluo- 
rescence photons (i.e., rp cc ). We use W = £f2 J Q T v(t)dt, 
which is £ times the work of the driving field per unit en- 
ergy defined by haj^, using the Bloch equation it is easy 
to show that when W >> 1 (i.e., photon count is large) 
W — ► £r J Q p cc {t)dt as expected. Using Eq. (||) the aver- 
age number of photons counted in time interval (0, T) is 
(n) = (W) where (• • •) is an average over the stochastic 
process. Mandcl's Q parameter is used to characterize 
the fluctuations and it is straightforward to show that 



(W 2 ) - (W)' 
(W) 



(4) 



We see that Q > indicating that photon statistics is 
super Poissonian. 

We now consider the important limit of weak laser in- 
tensity of Eq. (|2|) . We use a straightforward perturbation 
expansion in the Rabi frequency to find 



v = —Re 
2 



dii exp 



(•5) 



In standard line shape theories Eq. (||) is averaged, the 
well known result gives the Wiener-Khintchine formula 
for the line shape (7(ljl)) = hniT->oo(IU)/T, which is 
valid for stationary processes 



</K)) = i^-Re 



dTe -toLT-iv/a Ci (r) 



(6) 



where the one time dipole correlation function is Ci(r) = 
(e* f° ^3 A)IJ '' ^ dt ). Let us now consider the fluctuation, 

£2^4 ,T ,T ,T ,T 



dtidt 2 dt 3 dt4X 
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e -iu> h {t2-t 1 +t 3 -u)-r(\t 1 -t 2 \+\t 3 -u\)/2 ( j 3 / Ti)T2)T3 ^ ^ 

where the three time correlation function 
C 3 (7-1,7*2,73) = 



(exp 



i f 2 AwjtfW - 1 1 4 Y\ Awj(t')d*' 
Jtl Jt 3 



(8) 



contains the information on the stationary spectral diffu- 
sion process relevant for the calculation of Q. It depends 
on the time ordering of ^1,^2,^3,^4 for which there are 
4! = 24 options. In Eq. (||) t\ = tu - h, r 2 = tin - hi 
and 73 = <iv — tin and by definition ti < tu < tm < tiy. 
Let us consider the case t\ < t 2 < t 3 < ti (for this case 
t\ = ti,t2 = tu etc) and t\ = t 2 — t\, t 2 = t$ — t 2 and 
r 3 = i 4 — t 3 . We define the pulse function S(t) = — 1 for 
ti < t < t m S(t) = for t u < t < tm and S(t) = 1 
when tin < t < tjv, the shape of this pulse is shown 
in the first line of Table 1. Using this pulse the corre- 
lation function is written as a characteristic functional, 



C 3 (71,72,7-3) = (exp 



Simi- 



larly other pulse shapes describe the other time orderings 
as described in Table 1. We note that pulses 1 and 2 are 
encountered in theory of four wave mixing [fl3f . 

Eqs. (|4|-||) give a general prescription for the calcu- 
lation of Q for broad types of spectral diffusions. Now 
we investigate basic properties of line shape fluctuations 
considering a simple process |l2fl . We assume J = 1 and 
Au>(t) — vh(t) and h{t) is a random telegraph process 
hit) = 1 or h(t) = — 1. The transition rate from state 
up to down and vice versa is -R, generalizations for more 
complicated cases will be considered elsewhere. 

Using a method of Suarez and Silbey p0[ , developed 
in the context of photon echo experiments, we now ana- 
lyze the properties of the three time correlation function. 



Define the weig hts P? f (t) = (exp -is/c 'Au(t')dt' 



hf 
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where the initial (final) state of the stochastic process 
Au is i (/) and s = or s = ±1 or s = ±2 Q. For 

example P^_(t) is the value of (exp (i J Q * Auj(t')dt'^j) for 

a path restricted to have Acj(O) = v and Auj(t) = v. Us- 
ing these weights and based on the Markovian property 
of the process, we find for t\ < t 2 < t-j < 

where the summations are over all possible values of i — 
±, j = ± ,k = ±, and I = ±. 
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C 3 (Tl,T2,r 3 ) 


time order 
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|E4(n)^(r 2 ^(r 3 ) 


2314^4,2^3 



Table 1. Each pulse in column 1 represents four time 
orderings corresponding to a three time correlation func- 
tion in column 2. The symbol 1234i— 4,2— 3 corresponds 
to h < t 2 < t 3 < t 4 , U < t 2 < ta < h, h < i 3 < t 2 < t± 
and ti < t 3 < t 2 < fi. The complex conjugate (C.C) 
of C3(ti,T2,T3) in line 1 describes the time orderings 
21431-^4,2^3, C.C [C3 (ri, T2, r 3 )] of pulse 2 corresponds 
to 1243i^4,2^3, while C.C [C3 (n, r 2 , r 3 )] of pulse 3 cor- 
responds to 14231^4,2—3- 

In a similar way we consider the other 23 time order- 
ings. We find only three types of pulses S(t) are needed 
for the definition of the three time correlation function, 
these pulses are presented in Table 1 together with the 
corresponding C3(ti, t 2 , r 3 ). We then sum all the 24 con- 
tributions and use the convolution theorem of Laplace 
transform to arrive at the main result of this paper 

,2v r 1 r ^ V- 



(W 2 ) = £~ 1 { 



16m 2 



P^(u + u + )P°(u)P+ 1 (u + u_) 

+P ij 1 (u + u+ ) P% (u + r) P+ 1 (u + «_) 

+P+ 1 {u + u_)P%{u)P+ 1 (u + u^) 
+P+ 1 (u + U-)P% (u + r) P+ 1 (u + u-) 
+2P+ 1 (u + u-) P+ 2 (u + 2u_) P+ 1 (u + u-) + C.C}} 

(9) 

where u± = T/2 ± icjl, denotes the inverse Laplace 



transform and Pfj(u) is the Laplace transform of P"j(t). 
Using standard methods of complex analysis (one can 
also use standard numerical or symbolic packages) we 
find the exact analytical expression for Mandel's Q pa- 
rameter Eq. (0), based on Eqs. @ and (W) [i.e. Eq. 
(^) when T — ► 00]. It turns out that Q is not a simple 
function of the model parameters; however, as we show 
below in certain limits simple behaviors are found. 

We now consider the slow modulation limit i/>F> 
R. Using our exact result Eq. (Sh we find 



,-2RT 



2RT 



(A 



I-) 
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where I± = 2 r|E + i(uj^ =p z^) | 2 /4 are stationary solu- 
tions of time independent Bloch equation with frequency 
detuning (wl =p v) , respectively. Similar to the line shape 
(I (uj l )) ~ (J + +J_)/2, Q Eq. © exhibits splitting with 
two peaks at wl = ±za Unlike the line shape, Q depends 
on measurement time T and on R. As an example we 
choose the parameters F = 40MHz, ft = F/10, R = 1/scc 
and wl = v = 1GHz, which mimics a SM coupled to a 
slow two level system in glass ||, we find for RT ^> 1, 
Q = 2 x 10 5 £ which implies that fluctuations are very 
large (if compared with the fast modulation limit soon 
to be considered). The lcngthly though straightforward 
derivation of Eq. ([!(]) is based on careful analysis of the 
poles of (W 2 ) and (W) in Laplace domain. The result 
Eq. ( |To| ) can be easily understood in the following way. 
The molecule can be found in two states up + and down 
— , in the slow modulation limit the rate of photon emis- 
sion in these two states is determined by the stationary 
solution of time independent Bloch equation, namely I± , 
while transients are neglected. Using this physical pic- 
ture one can re-derive Eq. (|ic| ) in a straightforward way. 
Eq. ([h]) describes well the experiments from Bordeaux 
H of SMS in glass, there the second order correlation 
function g( 2 > (t) was used to characterize the dynamics of 
a SM coupled to a two level system in the glass. To see 
the relation between these experiments and our work we 
recall the relation fllq,19| 



Q = {iM) 



dt 2 g (2 \t 2 )-T 



and for slow process [t) = 1 + C exp(—2Rt) with C — 
(/+ — 1_) 2 /(/+ + /_ ) 2 , this expression being equivalent to 
Eq. (|Io|). The question now remains whether one can go 
beyond this slow modulation limit and investigate faster 
processes based on SMS. 

We therefore consider fast modulation limit for which 
only the T — > 00 limit is relevant. Mathematically differ- 
ent definitions of a fast modulation exist depending on 
the R — » 00 limiting procedure used. First consider the 
case R — > 00 keeping v and T fixed. This fast limit gives 
Q = 0. This is expected since in this case the bath is so 
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fast the molecule cannot respond to it (e.g., the line width 
of the molecule is T). A more physically interesting case 
is to let R — > oo but keep the ratio r e ff/r = v 2 / (RT) 
finite. In this fast modulation limit the well known line 
shape is Lorentzian with a width T + T c ff and as men- 
tioned it exhibits motional narrowing. Using our exact 
results we find in the same limit and for cJl = (i.e. at 
resonance) 



cff 



r(r+ r cS y (r + 2r eff ) 



(ii) 



To estimate the magnitude of these fluctuations we note 
that the maximum of Q is found when r cff = r(l + v / 5)/2 
and then Q ~ 0.36£fJ 2 /T 2 . Even if we take |fi|/r = 0.1 
and £=5xl0~ 2 asa reasonable estimate for weak 
laser field and detection efficiency, we find a small value 
Q ~ 1.8xl0~ 4 . However values of \Q\ < 5 xlO -4 were re- 
ported in |l6||2^ ] , so measurement in the fast modulation 
regime might also be possible. 

The fast modulation limit can be analyzed using a 
method developed by Loring and Mukamel Jl3| in the 
context of four wave mixing. Briefly, in this limit 
we use a factorization approximation Cs(tx,T2,T3) ~ 
Ci (n) Ci (r 2 ) Ci (r 3 ), and Ci (r) is calculated based on 
standard line shape theory (details will be given else- 
where p3f ). This approximation works well since corre- 
lation between the state of the molecule (i.e., + or — ) 
during one pulse interval [say S(t) = 1] with that of the 
following pulse interval [say S(t ) = 0] are unimportant, 
because the bath is fast. Employing this approach, or 
using Eq. (||) directly, we find 



Q 



gnT off 
2rK+r e V4) 



r<r, 



off 



(12) 



eg 



<r. 



Eq. ( |l2| ) exhibits several interesting behaviors (i) when 
r <C r c ff a simple relation between the line shape and 
Q holds, Q = 2(1 (o>l))/T, and both exhibit motional 
narrowing, (ii) when r e g <C T, Q ~ this behavior 
showing that as the bath becomes faster the deviations 
from Poisson statistics become smaller (similar behavior 
is observed in the slow modulation limit) (iii) a dip in 
Q is found when r c ff <C T (note that we have neglected 
0(r e g/r) 2 terms which give corrections at cjl = 0). 

To conclude, the three time correlation function de- 
scribes line shape fluctuations for a single molecule un- 
dergoing a stationary spectral diffusion process. We 
know that similar correlation functions yield important 
information in non-linear spectroscopy ]T^ | (i.e., beyond 
linear spectroscopy), similarly "Q SMS" yields detailed 
information on dynamical processes in the condensed 
phase (i.e., beyond the standard line shape theories). De- 
pending on the bath time scale, different time orderings of 



the three time correlation function contribute to the line 
shape fluctuations and in the fast modulation limit all 
time orderings contribute. Our exact results for a simple 
Kubo- Anderson process are applicable to the experimen- 
tal situation, in the slow and possibly in the intermedi- 
ate modulation regime (where Q is large). Our results in 
the more challenging fast modulation case give the the- 
oretical limitations on measurement accuracy needed to 
detect Q. 
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